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VHF SATELLITE NAVIGATION SYSTEM 


1.0 INTRODUCTION 

The ever-increasing need for an improved air traffic control 
system, particularly over the North AtLantic area, has sparred an 
intense search for system concepts and techniques which wiLL satisfy 
the future requirements of commercial air traffic This proposal 
describes a VHF Satellite Navigation Experiment which can be con- 
ducted using previously developed technology, equipments, and the 
ATS-1 and ATS- 3 satellites presently in orbit The primary goal of 
this experiment is to determine the position location accuracy that 
can be achieved and to define confidence levels that can be attained 
by ranging to geosynchronous satellites at VHF frequencies. A 
secondary goal is to demonstrate that aircraft positions can be deter- 
mined by these means in real time under flight conditions. 

Another potential application of this VHF Satellite Navigation 
System is location and tracking of balloon-borne meteorological 
platforms Recommendations and studies concerning the feasibility 
of a Global Weather Observing System indicate that an experiment 
to demonstrate such a position location capability would be mean- 
ingful and potentially important. This demonstration would be par- 
ticularly significant because the position location is accomplished via 
geosynchronous satellites. 

A limiting parameter which must be considered in the design 
of an air traffic control system, as well as a global meteorological 
balloon system, is the satellite bandwidth requirements per aircraft 
or balloon to be serviced As discussed in subsequent sections of 
this document, the proposed sidetone ranging technique requires a 
significantly lower satellite bandwidth than that required by other 
ranging schemes of comparable accuracy and resolution 
Additionally, the specific sidetone pattern proposed for this 
experiment represents a near-optimum configuration for sidetone 
ranging to geosynchronous satellites when considering ranging 
accuracy requirements, range ambiguity resolution, satellite 
bandwidth requirements, allowable acquisition time, etc Although 
sidetone ranging has been selected as the means for range measure- 
ments, the experiment has been designed to investigate fundamental 
limitations and statistical characteristics which are independent of 
the range measurement technique 

The VHF Satellite Navigation System described herein 
utilizes the technology and equipment developed for the Omega 
Position Location Equipment (OPLE) Experiment, and it employs 
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the Applications Technology Satellites (ATS-1 and ATS-3) for relay 
between the Ground Control Center (GCC) and cooperating aircraft. 

The GCC has the capability of determining the unambiguous Location 
of each cooperating aircraft, provided the aircraft is within the area 
of simultaneous coverage by the ATS-1 and ATS-3 satellites, and 
determination of the aircraft position is not dependent upon another 
system {such as LORAN or OMEGA) 

The data processing functions performed by the OPLE Control 
Center in the present OPLE Experiment are almost identical to the 
sidetone processing functions that are required of the GCC for the 
proposed experiment In. addition, the self-check tones presently 
generated by the OPLE Control Center represent a near-optimum 
tone configuration for sidetone ranging to synchronous satellites. As 
a result, the OPLE Control Center hardware and software can be used 
to perform this VHF Satellite Navigation Experiment with only a few 
modifications to the existing equipment being required. These modifi- 
cations can be made without affecting the functions that are required of 
the OPLE Control Center by the OPLE Experiment. 

This VHF sidetone ranging experiment is a logical extension of 
the OPLE Experiment Also, the sidetone ranging technique described 
herein is very similar to the ranging technique used in the proposed 
Position Location and Aircraft Communication Equipment (PLACE) 
Experiment described in Goddard Space Flight Center Document No 
X-731 -67- 159, thus, the results obtained in this VHF sidetone ranging 
experiment can be incorporated into the design of the PLACE system. 
Just as with the existing OPLE system hardware, this experimental 
VHF Satellite Navigation System is not designed to be directly 
expandable into an operational system, however, the information pro- 
vided by this experiment will most certainly be useful in. the subse- 
quent design of any operational air traffic control system 

Position location m the proposed system is accomplished by 
means of simultaneous ranging from two geosynchronous satellites 
to an aircraft of known altitude. This ranging scheme utilizes the 
unmodified VHF sidetone pattern that is presently generated m the 
OPLE Control Center for self-check purposes From the measured 
ranges between the aircraft and the ATS-1 and ATS-3 satellites, 
two spheres are defined. These spheres are centered at the two 
satellite positions and they include the aircraft position Similarly, 
a sphere is defined from the known aircraft altitude relative to the 
' center of the Earth. This third sphere is concentric with the Earth 
and also includes the aircraft position. The intersections of these 
three spheres define two possible aircraft positions which are 
symmetrical about the Earth's equator. The position ambiguity 
between the Northern and Southern hemispheres can be readily 
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resolved from a priori knowledge for all flight paths other than those 
near or crossing the Earth’s equator. For these cases, position 
ambiguity can be resolved by continuous tracking of the aircraft 
position from a known starting point 

Range measurements from the aircraft to the two ATS 
satellites are performed simultaneously in this experiment; however, 
only one sidetone pattern (one data channel) is transmitted to the 
ATS-1 satellite during a range measurement, This is mandatory to 
avoid the inter modulation distortion which results from the hard 
limiter in the ATS-1 satellite transponder when a multiplicity of sig- 
nals exists. As a result, the aircraft terminal can be channelized to 
receive on a given channel frequency and to transmit at another chan- 
nel frequency. In an operational system, this would minimize the 
bandwidth requirements for each aircraft, thereby maximizing the 
number of aircraft that could be serviced simultaneously with a par- 
ticular satellite bandwidth. Additionally, only one transponder is 
required in the aircraft terminal to perform simultaneous range 
measurements to two satellites, thus, the terminal hardware require- 
ments are kept to an absolute minimum This hardware reduction' 
would be of utmost importance for data collection platforms to be used 
with meteorological balloons. 

This experimental VHF ranging system consists of (1) a modi- 
fied OPEE Control Center, including a second VHF antenna, which 
functions as the Ground Control Center, (2) an aircraft terminal, and 
(3) a reference terminal and satellite transponder simulator which 
will be used for system calibration purposes. The system will also 
utilize the VHF transponders m the ATS-1 and ATS- 3 satellites. The 
GCC will be capable of receiving and processing either VHF ranging 
data or OPLE data, therefore, the present OPEE Experiment and the 
proposed VHF Satellite Navigation Experiment can be conducted con- 
currently if desired; however, a brief down time is required to switch 
from one experiment to the other 
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PRECEDING PAGE BLANK NOT FILMED. 

2. 0 EXPERIMENT DESCRIPTION 

The primary objective of this VHF Satellite Navigation 
Experiment is to determine the accuracy and reliability Limitations 
of Locating the positions of aircraft (or remote platforms) by means 
of VHF sidetone ranging measurements to synchronous sateLLites 
The range measurements are made from a Ground Control Center 
(GCC) Located at Goddard Space Flight Center to cooperating air- 
craft via VHF links through the ATS-1 and ATS-3 sateLLites. The 
resulting range data is combined with known aircraft attitude 
information to Locate the position of the aircraft. The experiment 
can be performed with an FAA aircraft, commercial air freight 
aircraft, or any other suitable aircraft that is available. Aircraft 
performance characteristics, such as speed, are not limiting 
factors. 


The experiment configuration is shown m Figure 2-1. The 
GCC transmits VHF ranging tones to the ATS-3 satellite which 
relays the tones back to the Ground Control Center and to cooperating 
aircraft. The aircraft terminal receives the ranging tones from the 
ATS-3 satellite and retransmits them simultaneously to the GCC via 
the ATS-1 and ATS-3 satellites, The GCC receives and demodulates 
the ranging tones from the ATS-1 and ATS-3 satellites and compares 
the delay of each to the delay of the ranging tones that were trans- 
mitted by the GCC The delay differences between the received and 
transmitted tones represent the various VHF path lengths between 
the GCC, the ATS satellites, and the aircraft. These simultaneous 
range measurements, along with a similar but separate range 
measurement from the GCC to the ATS-1 satellite, are used to 
determine the ranges from the aircraft to the ATS-1 and ATS-3 
satellites. 

As indicated above, this experiment employs the VHF trans- 
ponders in the ATS-1 and ATS-3 satellites. Each transponder has 
transmit and receive frequencies of 135 6 MHz and 149.22 MHz, 
respectively Similarly, a VHF transponder is required m the air- 
craft terminal, although it is inverted to provide transmit and 
receive frequencies of 149-22 MHz and 135 6 MHz, respectively. 

The GCC consists of the present OPLE Control Center, slightly 
modified for this experiment, however, the functions required of 
the OPLE Control Center by the OPLE Experiment are retained 
This permits concurrent performance of the VHF Satellite Naviga- 
tion Experiment and the OPLE Experiment, if desired. 

The ranging method proposed for this VHF Satellite Naviga- 
tion Experiment is the sidetone ranging technique. Three non- 
harmonic tones of appropriate frequency are transmitted with a 


5 




Figure 2-1, VHF Satellite Navigation System Configuration 
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reference carrier and then used m such a manner that the desired 
ranges can be determined unambiguously. The ranges measured 
between the aircraft and the ATS-1 and ATS-3 satellites define two 
spheres centered at the satellite positions, and a third sphere is 
defined by the known aircraft altitude above the center of the Earth. 

The position of the aircraft is determined by means of digital computer 
processing of the measured VHF range data and the known aircraft 
altitude information. 

Another possible application of the VHF Satellite Navigation 
System involves meteorological balloons that are currently being 
launched from Palestine, Texas, by the National Center for Atmo- 
spheric Research (NCAR), A balloon-borne terminal could be launched 
in the January-March 1969 time period when prevailing westerly 
winds would place the balloon in simultaneous view of the ATS-1 and 
ATS-3 satellites, and the balloon could then be tracked by the GCC to 
measure the wind velocity at the balloon's flight altitude (approximately 
20 kilometers) . 

The absolute position location accuracy of this VHF sidetone 
ranging system is determined primarily by the uncertainty of the iono- 
sphere. Improved ionospheric prediction techniques are currently being 
developed, however, by using presently available ionospheric data to 
compute bias errors, it is estimated that a worst-case absolute posi- 
tion location accuracy on the order of ±2. 0 kilometers (±1. 08 nautical 
miles) can be achieved. For a differential (relative) position location 
measurement, such as that required in an air traffic control system, 
the ionospheric bias errors are common. In this case, the estimated 
worst-case error is on the order of ±0. 7 kilometer (±0.3 8 nautical 
mile), provided the distance between- the two aircraft is withm approxi- 
mately 100 kilometers. 
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3. 0 SYSTEM DESCRIPTION 

Position location of an aircraft (or remote platform) in the 
YHF Satellite Navigation System is accomplished by first determining 
the distance from the aircraft to the center of the Earth and the ranges 
from the aircraft to the ATS-1 and ATS- 3 satellites This altitude 
and range data is then processed off-line to (1) resolve the range 
ambiguities, (2) correct for system bias errors, and (3) 'determine the 
aircraft position. Controlled portions of this experiment will be per- 
formed with the aircraft altitude known to good precision. Hence, the 
distance from the aircraft to the Earth’s center can be determined 
by an iterative process as described later, and other factors contrib- 
uting to position uncertainties can be investigated separately. Therefore, 
the primary function of the YHF Satellite Navigation System in this 
experiment is to measure the aircraft-to- satellite ranges and to 
compute the aircraft position 

The distances from the aircraft to the ATS satellites are 
measured indirectly by means of a YHF sidetone ranging technique. 

In this system, a set of three non-harmonic VHF ranging tones and 
a reference carrier are generated at the GCC and transmitted to the 
ATS-3 satellite at 149. 22 MHz. The ATS- 3 satellite translates this 
ranging tone pattern to 135. 6 MHz, maintaining the original tone 
phase relationship, and relays the ranging tones back to the GCC and 
also to an aircraft The return ranging tones are received and 
processed by the GCC to determine the unambiguous two-way range 
from the GCC to the ATS-3 satellite; this range is represented by 
the differences in phases of the received and transmitted ranging 
tones at the GCC. The aircraft also receives the 135. 6-MHz ranging 
tone pattern from the ATS-3 satellite and translates it to 149. 22 MHz 
without altering the tone phase relationships. The ranging tones are 
then relayed back to the GCC via the ATS-1 and ATS-3 satellites 
The GCC receives and processes the two tone patterns from the 
satellites to determine (1) the unambiguous two-way range from the 
GCC to the aircraft via the ATS-3 satellite, and (2) the unambiguous 
one-way range from the GCC to the aircraft via the ATS-3 satellite 
and the return from the aircraft to the GCC via the ATS-1 satellite. 

In a similar but separate measurement, the GCC transmits the 
sidetone pattern to the ATS-1 satellite. The return signal from the 
satellite is then received and processed by the GCC to determine 
the unambiguous two-way range from the GCC to the ATS-1 satellite. 
Finally, the distances from the aircraft to the ATS-1 and ATS-3 
satellites are computed using range differences derived from the 
above series of range measurements. 

The computed distances from the aircraft to the ATS satellites 
are processed off-line to remove the ranging errors due to 
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ionospheric effects and other system bias components. These corrected 
ranges are then combined with known aircraft aLtitude data to obtain the 
aircraft position. 

This VHF Satellite Navigation Experiment employs a Ground 
Control Center, a reference terminal, an aircraft terminal, the ATS-1 
and ATS-3 satellites, a satellite transponder simulator, an off-line 
digital computer, and at least one cooperating aircraft. The present 
OPLE Control Center can be modified and used as the GCC, whereas 
the aircraft and reference terminals can be constructed from modified 
commercial equipment. The remaining equipment will be available fox 
use in the experiment. 


3 . 1 Ground Control Center (GCC) 

The GCC must generate the required VHF ranging tones, 
transmit these tones to the ATS-3 satellite or ATS-1 satellite as 
required, and then receive and process the return signals to determine 
the phase differences between the transmitted and received ranging 
tones. In addition, the GCC must perform all the functions that are 
presently required by the OPLE Experiment. It will be shown that- the 
additional capability can be provided by relatively minor modifications 
of the existing OPLE Control Center (OCC). 

3.1.1 Ranging Tone Generation 

The self-check subsystem in the GCC generates a channel of 
simulated OPLE Platform Electronics Package (PEP) data. This simu- 
lated PEP data includes an acquisition/reference tone (carrier) and 
three frequency-translated Omega tones (Figure 3-1) which function 
as the ranging tones in this VHF ranging experiment. In the OPLE 
Experiment, these tones are commutated in synchronization with the 
Omega Navigation System timing and then transmitted to the ATS-3 
satellite; the return signal from the satellite is received and processed 
by the GCC to determine the phase differences between the received and 
transmitted tones. This is identical to the tone processing that is 
required in this proposed VHF Satellite Navigation System. 

To enable VHF sidetone ranging measurements via the two 
satellites, the GCC self-check subsystem must be modified to provide 
an alternate mode of operation in which the three simulated Omega tones 
are continuous rather than commutated. This modification consists of 
relatively minor wiring changes. 

3.1.2 RF Electronics Subsystem Modification 

As described in Section 3. 0, continuous range measurements 
to an aircraft are performed by transmitting ranging tones to the 
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Side tone Ranging Pattern 


ATS-3 satellite -while simultaneously receiving return ranging tones 
from the ATS-1 and ATS-3 satellites. Hence, dual VHF antennas and 
receivers must be incorporated into the GCC as shown in Figure 3-2. 
The return ranging tones that are received by direct translation from 
the ATS-3 satellite are assigned a channel frequency different from 
that used to relay the ranging tones to the aircraft. However, the two 
sets of ranging tones received at the GCC from the aircraft via the 
ATS-1 and ATS-3 satellites occupy the same channel, with the neces- 
sary isolation between these two signals being provided by the GCC 
antennas. Additional isolation is provided because of the difference 
in the two satellite transponder "side-step" frequencies. Because of 
this relatively large difference in the “side-step" frequencies of the 
two satellites, separate VHF frequency synthesizers must also be 
included for the two VHF receivers in the GCC. 

The GCC receivers extract the relative phase data from the 
return ranging tones. The first step of this extraction process is to 
remove the random phase shifts which result from the series string 
of frequency translation elements m the VHF paths. These random 
phase shifts apply equally to all ranging tones, including the carrier, 
thus their effect can be removed by phase-locking to the reference 
tone. This is accomplished in the OP1.E receivers that follow the 
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VHF receivers. .After the phase jitter has been removed, the phases of 
the three received ranging Tories are compared to the phases of the 
ranging tones transmitted by the GCC, and the resulting phase angles 
constitute the basic ranging data. 

For the separate range measurement to the ATS-1 satellite, 
the ATS-3 antenna at the GCC is positioned to point at the ATS-1 
satellite and the ranging tones are transmitted to it. The return tones 
from the ATS-1 satellite are received at the GCC by the ATS-1 antenna 
and receiver because the desired satellite "side-step" frequency error 
compensation has already been introduced into this receiver. 

A functional GCC block diagram illustrating the required 
OPLiE Control Center modifications is shown in Figure 3-2, with the 
added equipment enclosed in the dashed lines. It is significant to note 
that the ATS-3 receiver at the GCC requires a pilot tone tracking 
phase-lock loop to remove the excessive frequency jitter that is 
inserted by the ATS-3 satellite. This is not necessary in the ATS-1 
receiver at the GCC, as was demonstrated during the successful OPLE 
Engineering Tests with the ATS-1 satellite. 

3.1.3 Timing and Control Modifications 

Because the system timing and data processing requirements 
are slightly different for the VHF ranging mode of operation, minor 
modifications to the control electronics subsystem in the GCC will be 
required. These modifications consist only of wiring changes. The 
OPLE system software will not have to be changed for this experiment, 
however, a modification to it would allow a reduction m the ranging tone 
processing time, 

3.1.4 Off-Line Data Processing 

The off-line data reduction and data processing required to 
locate an aircraft by means of VHF sidetone ranging is entirely 
different from that required for position location by means of the OPLE 
technique. Consequently, additional off-line software must be developed 
for the VHF Satellite Navigation Experiment. 

3.1.5 OPLE Experiment Compatibility 

The GCC described herein is an OPLE Control Center with 
the tiecessary modifications included to provide the additional capa- 
bility required to perform the VHF ranging experiment along with the 
OPLE Experiment. Equipment, wiring, and software changes are 
minor, and they can be made in such a way that only a few minutes 
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down time will be required to switch from one experiment to the 
other. The GCC can be used to perform both experiments, thus pro- 
viding the capability of making nearly simultaneous position location 
measurements for comparison of the two position determining techniques. 


3 . 2 Reference Terminal 

The reference terminal is comprised of a YHF transponder 
and a YHF antenna, it can be used by the GCC operator to verify proper 
operation of the VHF SatelLite Navigation System and to perform system 
calibration without relying on the aircraft terminal. The reference 
terminal can be located near the GCC or it can be installed at a remote 
ground facility; however, in either case, adequate coverage by both 
ATS satellites must be provided. This terminal can be used at the 
remote location to perform differential position location measurements 
in conjunction with the aircraft terminal, and it can also be used for 
ionospheric error studies. 

A block diagram of the reference terminal transponder is 
shown m Figure 3-3. It consists of a double - con ve rs ion VHF receiver 
which receives and amplifies the VHF ranging tones and pilot tone 
that are relayed from the GCC via the ATS-3 satellite. The ranging 
tones and pilot tone are separated by band- pass filters. Then, a two- 
step phase jitter suppression operation follows to completely remove the 
phase jitter that is introduced in the ranging tones by the ATS-3 satellite 
and the VHF links. The necessity for doing this has already been demon- 
strated in the OFLE Experiment, and the required OPEE Control Centex 
modifications have been made to accommodate the frequency jitter 
acquired in a one-way pass through the ATS-3 satellite. In the first 
jitter suppression stage, phase -lock to the pilot tone is acquired by a 
wideband phase -lock loop, thus the group phase j'itter is tracked out. 

This includes any doppler shift frequency, satellite "side - step" fre- 
quency error, satellite frequency jitter, or other phase jitter. The 
ranging tones are next applied to a narrowband phase -lock loop which 
acquires lock to the reference tone and thus tracks out any phase jitter 
on the reference tone relative to the pilot tone. The jitter -free ranging 
tones are then up- converted and transmitted at a power level of 100 
watts. To prevent saturation of the transponder receiver by the power 
amplifier, approximately 70 dB of isolation is provided by the pre- 
selector between the antenna and the receiver. Because the VHF trans- 
ponder may be required to operate in the vicinity of high-power 
transmitters, a high level of image rejection is necessary, hence the 
double -conversion receiver technique is employed. 
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Figure 3-3. BLock Diagram of Reference Terminal Transponder 
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The reference terminal will 
be constructed from commercial 
equipment with the exception of the 
phase-lock loops, which will be of 
the type presently used in the OPLE 
Control Center. 


3 . 3 Aircraft Terminal 

The aircraft terminal for 
the VHF Satellite Navigation Sys- 
tem. consists of a VHF transponder 
(which is a ruggedized version of 
that used in the reference termi- 
nal) and a VHF aircraft antenna 
(Type DMC 33-2). This equipment, 
shown in block diagram form m 
Figure 3-4, can be installed in a 
cooperating aircraft. The opera- 
tion of the transponder is as out- 
lined in Section 3.2. 

The aircraft can also be 
equipped to operate m the OPLE 
mode if desired. In this OPLE 
mode of operation (VHF sidetone 
ranging not performed), a PEP is 
utilized as the aircraft terminal. This mode of operation is identical to 
the OPLE system operation m that the PEP (1) responds to coded inter- 
rogations from the GCC via the ATS-3 satellite, (2) encodes and trans- 
mits sensor data to the GCC via a return link through the satellite, 

(3) receives and conditions the VLF tones transmitted by the Omega 
Navigation System, and (4) up-converts the Omega tones to the VHF 
link frequency and transmits them to the GCC via the ATS-3 satellite. 


VHF 

ANTENNA 



572 3 7 


Figure 3-4. Aircraft Terminal 
Configuration 
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4.0 


VHF SIDETONE RANGING 


Location of the position of an aircraft in this navigation system 
depends on four VHF range measurements. The range measurements 
are shown in Figure 2-1 and include: 

(1) The path length from the GCC to the ATS-3 satellite 
and back (2 Rgj) 

(2) The path length from the GCC to the aircraft and back 
via the ATS-3 satellite (2Rg3 a ) 

(3) The path length from the GCC to the aircraft via the 

ATS-3 satellite and back to the GCC via the ATS-1 

satellite and back (R„ Q „ i „) 

gialg' 

(4) The path length from the GCC to the ATS-1 satellite 
and back (2Rgj). 

The first three range measurements are made simultaneously and the 
fourth is made on a time- shared basis. The primary justification for 
time-sharing the range measurements is that potential multiple access 
difficulties due to the ATS-1 satellite hard-limiting transponder can be 
avoided. Furthermore, the ATS-1 satellite is nominally stationary and 
its dynamic behaviour is such that frequent measurements of the range 
from the GCC to the ATS-1 satellite are completely adequate. This 
approach also allows the use of minimal hardware m the aircraft trans- 
ponder and at the GCC. 

Since the range measurements are made at VHF, the range 
estimates are perturbed by ionospheric effects. However, it is possible 
to correct for some of these errors, thus making a VHF Satellite Navi- 
gation Experiment possible with relatively high accuracy being achieved 
A functional block diagram showing the operations leading to a VHF 
range measurement is shown in Figure 4-1. 

The ranging method used in this system is the sidetone ranging 
technique. This is one of the two techniques usually considered for 
ranging between a fixed station and a cooperative target; the other 
standard ranging technique employs a pseudo-noise (PN) code. In 
Appendix A, the relative merits of these two ranging techniques are 
discussed using the criteria proposed by Baghdady and Kruse (Refer- 
ence 1). It is shown conclusively that sidetone ranging is the preferred 
technique for this application. 

4. 1 VHF Range Resolution 

The range value for each of the four measurement paths in this 
VHF Satellite Navigational System is determined from simultaneous 
phase measurements of the three sidetones (426 Hz, 707Hz, and 941 Hz). 
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Figure 4-1. Functional Block Diagram, of VHF Range Measurement 






These phase measurements represent the phase shifts of the respective 
side tones due to their transmission through the path of unknown range. 
The range (R) is related to the phase shift (0) and the sidetone frequency 
(f) by the equation: 


R 


c0 

£ 


where c is the propagation velocity. 


The range resolution is limited by the phase measurement 
accuracy and is inversely proportional to the sidetone frequency. A plot 
of the RMS range error due to noise for each of the sidetone frequencies 
is shown in Figure 4-2, Included in this figure is a plot of the range error 
which results from the best weighted-average of the three independent 
range measurements. A typical value of phase signal- to-noise ratio for 
this system is 37 dB, yielding RMS range errors of 520, 320, 230, and 
180 meters for the 426-Hz tone, 707-Hz tone, 941-Hz tone, and the best 
weighted-average, respectively. Factors which contribute to the phase 
noise include atmospheric noise, multipath effects, and system calibra- 
tion inaccuracies. A detailed discussion of these various factors, with 
their resulting range errors, is presented in Section 6. 1. 


The best weighted-average of the three measurements is com- 
puted to minimize the range error. This minimum error is attained by 
applying weighting functions proportional to the measurement accuracies, 
which in the case of constant phase jitter due to noise, becomes propor- 
tional to the sidetone frequencies. In other words, the gene ral weight- 
averaged range is expressed by the equation: 


R 



K, 


#2 

4-K 2 + 



Ki + K 2 + K 3 


% 


where the optimum vaLues of K are: 


K 1 

K- 


K, 



With these weighting vaLues, 

R = 


the equation becomes: 


c 


£ i + h + f 3 


This best weighted-average gives an equivalent 2 dB improvement over 
the accuracy that could he attained from the 941-Hz measurement alone. 
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RANGE MEASUREMENT ERROR (KM) 



Figure 4-2. RMS Range Error Due to Phase Error as a Function of SNR 
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4.2 


Range Ambiguity 

For each phase measurement made at the GCC on the returned 
sidetones, there exists a range ambiguity due to an unknown number 
of integral cycles associated with the total path length. Therefore, each 
of the range measurements made at the three sidetone frequencies has 
a systematic ambiguity. This range ambiguity can be resolved in this 
sidetone ranging system by off-line digital processing. This involves 
combining multiples of the measured phase values to produce a set of 
calculated phase values which correspond to synthesized sidetones at 
lower frequencies. A sufficiently low sidetone frequency is used such 
that the range can be defined unambiguously within the constraints of 
satellite /aircraft geometry. Once the unambiguous phase value is cal- 
culated, a "walk-up" method is used to sequentially resolve the 
ambiguity of each phase estimate, or measurement, at the next higher 
sidetone frequency. 

The frequencies of the sidetones generated at the GCC are 
426.4 Hz, 706.9 Hz, and 941.4 Hz, and a phase measurement (ambi- 
guous in the number of integral cycles) is made at each sidetone fre- 
quency. These phase measurements can be combined by means- of 
off-line processing to produce calculated phase values at other 
frequencies because the frequencies of the generated sidetones are 
non -harmonically related. One such set of synthesized frequencies 
(which will be used m this experiment) is calculated from the follow- 
ing combinations of the measured sidetone phase values: 


234. 5 Hz 
145. 9 Hz 
88. 6 Hz 
46. 0 Hz 
31.3 Hz 
14.7 Hz 


[941.4 - 706. 9] Hz ' 

[2(426.4) - 706. 9] Hz 
[941.4 - 2(426.4)] Hz 
[2(706.9) - 941.4 - 426.4] Hz 
[2(941.4) + 706.9 - 6(426.4)] Hz 
[706. 9 - 3(941.4) + 5(426.4)] Hz . 


The lowest synthesized tone of 14. 7 Hz can be used to define 
a range unambiguously up to a path length of 20, 000 Km, which is more 
than sufficient for resolving the unknown portion of the range over the 
various measurement paths in this experiment. The maximum unknown 
range value (approximately twice the Earth's radius, or 12,756 Km)' 
occurs in the two-way path between the ATS-3 satellite and the aircraft 
terminal. 
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Because the range resolution capability of the 14. 7-Hz tone 
is limited, a "walk-up” method is applied to sequentially resolve the 
ambiguity in the synthesized tones while proceeding to the ultimate 
range resolution provided by the three basic tones. The probability of 
incorrect ambiguity resolution as a function of signal-to-noise ratio of 
the generated tones is shown in Figure 4-3 for the three lowest synthe- 
sized tone values of 46. 0 Hz, 31.3 Hz, and 14.7 Hz. The significance 
of the three lower -frequency tones is that the 3 1.3 -Hz tone can he used 
to resolve range ambiguity over 84 percent of the possible path range 
with 5 dB less signal amplitude than the 14. 7-Hz tone for equal perfor- 
mance. Similarly, the 46. O-Hz tone can be used to resolve ambiguity 
over 57 percent of the range and thus reduces the signal amplitude 
requirement another 4.5 dB. Other methods of resolving the total ambi- 
guity would have to be applied with either of these higher-frequency tones 
however, in many cases this would not be difficult because the 46. O-Hz 
tone would provide a minimum value of range ambituity of 5, 000 Km 
along the Barth's surface. 

Some improvement in the probability of correctly resolving 
range ambiguities could be realized by adjusting the relative ampli- 
tudes of the three transmitted tones; however, this is not propo'sed in 
this experiment because it would require additional modification to the 
OPLE Control Center. 

4 . 3 VHF Link Considerations 

The VHF links that must be considered in the power budget 
analysis for this VHF ranging experiment are listed below m order of 
increasing signal-to-noi$e ratio; 

(1) The link from the GCC to the aircraft and return via the 
ATS-3 satellite 

(2) The link from the GCC to the aircraft via the ATS-3 
satellite and return via the ATS-1 satellite 

(3) The link from the GCC to the ATS-3 satellite and return 

(4) The link from the GCC to the ATS-1 satellite and return. ^ 

This order is dependent upon the fact that both the ATS-1 and the ATS-3 
are in full view of the aircraft and the GCC. In fact, the ATS-1 suffers 
an additional signal attenuation of the GCC due to its low elevation angle 
(3°); this degradation is approximately 5 dB. However, a higher-gam 
VHF antenna (18 dB rather than 13 dB) will be used for the ATS-1 
receiver at the GCC. The return from the aircraft via the ATS-3 
satellite has the least margin because the ATS-3 satellite transponder 
must relay a greater number of signals than the ATS-1 satellite 
transponder. 
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Figure 4-3. Ambiguity Resolution Versus Phase Noise 
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In the link analysis presented in the following paragraphs, 
optimum power sharing of the ATS-3 satellite power is first computed 
to maximize the performance of link (1). The performance of the other 
three links is then computed based upon these parameters. The worst- 
case analysis is then made without any power readjustment. 

t 

4.3.1 Optimum Ratio of GCC Power to Aircraft Transmitter Power 

In this section, the ATS-3 satellite power sharing is determined 
so as to minimize the total noise power density that is added to the rang- 
ing tones when they are relayed from the GCC to the aircraft and back 
via the ATS-3 satellite. This optimization is accomplished by noting the 
following facts: 


(1) 


The satellite is power limited (P s watts), and this power 
must be shared between (a) the tones relayed from the 
GCC to the aircraft (P sa watts m 4 tones), (b) the 
pilot tone relayed from the GCC to the aircraft (Pj 
watts); (c) the ranging tones and aircraft receiver 
noise relayed from the aircraft to the GCC (P g g + 
P S g n watts); and (d) the satellite noise that ismot sup- 
pressed due to the sum of signals a, b, and c (P w 
watts). It will be shown that P w is negligible, and 
hence the power summation is: 


P 


s 


P sa + 


P + P + 
sg sgn 



The reradiated aircraft receiver noise (Pgg^j is due to 
the aircraft transponder not being completely captured 
by the received ranging tones from the ATS-3 satellite. 
This power term can be kept arbitrarily small by a 
more complex terminal transponder design. For the 
conventional single-lF tone filter proposed for the 
transponder, however, this term may not be negli- 
gible. The ratio of the noise term (P S g n ) to the signal 
term (Psg) I s determined by the signal-to-noise ratio 
in the aircraft terminal transponder and hence, 

N a E a ' 

P sa L as 

where N a = noise power density in the terminal 
transponder receiver 

B a = terminal transponder tone filter bandwidth 
L as = link loss between aircraft and satellite. 


P sg + P sgn “ P sg 
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Nominal 

Worst Case 

Units 

Satellite Transmitter 




Coupling Loss 

-1.8 


dB 

Antenna Gam 

8. 5 


dB 

Path Loss . 

-167. 0 

-168. 1 

dB 

Polarization Loss 

-3 

-4 

dB 

Multipath Loss 

0 

-4 

dB 

Aircraft Receiver 

Antenna Gain 

2 

. 0 

dB 

Coupling Loss 

-2 


dB 

Link Loss 

-163.3 

-171. 4 

dB 

Antenna Temperature 

1000 

1600 

°K 

Receiver Noise Figure 

3 

4 

dB 

Effective Receiver 

1300 


° K 

Temperature 

Receiver Noise Power Density 

-197. 5 

-195. 5 

dBw /Hz 

Receiver Noise Power Density 


-34. 2 

-24. 1 

dBw /Hz 


3. 8 X 10~' 

i 3.9X10" 2 3 

watts /Hz 


Figure 4-4. Noise Power Density/Link Loss (N a /L as ) — Satellite to Aircraft 

(Frequency = 13 5.6 MHz) 


Using the nominal value for N a /L as (from Figure 4-4) 
of 3. 8 X 10" 4 watts /Hz, an IF bandwidth of 4 kHz, and 
an estimated P sa of approximately 10 watts, the nominal 
value of Psgn i- s negligible. Hence, the satellite power 
summation for nominal conditions becomes; 


P J = 


P + P 
sa sg 


(2) The contribution to the total noise by the GCC-to-satellite 

and aircraft-to-satellite links can be kept arbitrarily 
small by increasing the GCC and aircraft transmitter 

powers. Hence, the ratios of total noise power density 
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to signal power for the four ranging tones received at 
the 'GCC can be written as: 


N 

P 


N, 


P sg L gs 




^sa^as 


where Ng = noise power density- in the GCC receiver, 

Lg S = link loss between the GCC and the satellite. 

The other terms have been previously defined. Substitut- 
ing for Pgg from item (1), differentiating with respect to 
P ga , and solving the resulting quadratic equation for 
that value of P sa giving the minimum value of N/P, 
gives: 


1 - 


sa 


= (P fi - 


p j> 


fN L 
g a- s 

1 


1 /2 T 


N L 

g as 

L N 
sg a 


It has been observed in OPLE tests with the ATS-3 
satellite that the pilot tone must acquire approximately 
25 percent of the total satellite power to achieve 
reliable operation. For P g = 40 watts and using the 
nominal N/P ratios calculated m Figures 4-4 and 4-5, 
the optimum value of P sa is: 

P -22 watts 
sa 

and P„„ = 8 watts, 

sg 

(3) Because the relative GCC and aircraft transmitter power 
levels control the satellite power sharing, the following 
equation can be written, subject to the assumption made 
in items (1) and (2): 


or 


, L g-g P g_ 


sa 


L as P a 


sg 


P g _ ^as P sa 

P ~ P L 
a sg gs 
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Nominal Worst Case Units 


Satellite Transmitter 


Coupling Loss 

-1.8 


dB 

Antenna Gain 

8. 5 


dB 

Path Loss 

-167.0 

-168. 1 

dB 

Polarization Loss 

-3 

-4 

dB 

Ground Receiver 

Antenna Gain 

13 


dB 

Coupling Loss 

-2. 5 


dB 

Link Loss 

-152, 8 

-154.9 

dB 

Antenna Temperature 

1300 

1600 

°K 

Receiver Noise Figure 

3 _ 

4 

dB 

Effective Receiver 

1600 

2050 

° K 

Tempe rature 

Receiver Noise Power Density 

-196.6 

-195. 5 

dBw /Hz 

Receiver Noise Power Density /Link 

-43.8 

-40. 6 

dBw/Hz 

Loss 


4. 17 X 10 

-5 8.7X10 -5 

watts /Hz 


Figure 4-5. Noise Power Density/Link Loss (Ng /L S g) —Satellite to Ground 

(Frequency = 135.6 MHz) 

where = GCC transmitter power in the tones 

P a = aircraft transponder power. 

Again, substituting nominal values: 



Hence, if P a is chosen to be 100 watts, Pg is 25 watts. 
This is well within the capabilities of the present OCC 
transmitter . 
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4.3.2 Link Power Budgets 


In the previous section, the optimum satellite power division 
and GCC and aircraft transmitter power levels were calculated by making 
a number of simplifying assumptions. In this section, the actual nominal 
and worst-case ratios of noise power density to signal power per ranging 
tone are calculated for the various ranging paths, using information 
presented in Reference 7. First, the following accurate results are 
required. Using the derivation in item (2) of the last section, but not 
making the simplifying assumptions used in that section: 


P 


sa 


Ps 
1 + 


jlLl 


L g= p 


g 


Using the exact result for item (3) in the last section: 

P t 


T P P 
as a sa 


sg 


L gs P g 


1 + 


N B 
a a 

■^sa^as 


The nominal values, using B a = 4 kHz and the results of Figures 4-4 and 
4-5, are; 


P ga = 22 watts 
P S g = 8 watts. 

The worst-case values, using the results in the same figures, are: 

P = 27. 6 watts 
sa 

P S g = 1 . 56 watts. 

4.3.2. 1 Link from GCC to Aircraft and Back Via ATS-3 Satellite 

The effective noise power ratio per ranging tone is given by: 





/ \ 


’n" 


N N 

f N B 1 

i N N 

- 4 

— + s 

1 + a a 

U a + 3 

_p_ 

1 

P sg L g S ^a-^as | 

( P saL aS j 

* p sa p as ^g^gs 


Inserting nominal values from Figures 4-4 through 4-7 gives: 



1. 2 X 10“ 4 per Hz, or -3 9. 2 dB/Hz. 
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Nominal 


Worst Case 


Units 


Aircraft Transmitter 


Coupling Loss 

-2 


dB 

Antenna Gain 

2 

0 

dB 

Path Loss 

-168. 5 

-169. 7 

dB 

Polarization Loss 

-3 

-4 

dB 

Multipath Loss 

0 

-4 

dB 

Satellite Receiver 




Antenna Gain 

8 


dB 

Coupling Losses 

-1.8 


dB 

Link Loss 

-165.3 

-173. 5 

dB 

Antenna Temperature 

800. 


° K 

Receiver Noise Temperature 

400 

1200 

° K 

Effective Receiver 

1200 

2000 

° K 

Tempe rature 


• 


Receiver Noise Power Density 

-197. 8 

-195. 5 

dBw/Hz 

Receiver Noise Power Density/Link 

-32. 5 

-22 

dBw /Hz 

Loss 

5. 6 X 10' 

‘ 4 6 .3 X 10“ 3 

watts /Hz 


Figure 4-6. Noise Power Density /Link Loss (N s /L ag ) — Aircraft to Satellite 

(Frequency = 149, 2 MHz) 

Using simultaneous worst-case values, however: 

N = 1 . 2 X 10“ 3 per Hz, or -29. 2 dB/Hz. 

LpJi 

4. 3. 2. 2 Link from GCC to Aircraft Via ATS-3 Satellite and Return 
Via ATS-1 Satellite 

This link will have slightly less margin than the link described 
in Paragraph 4.3. 2. 1 
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Nominal 

Worst Case 

Units 

GCC Transmitter 




Coupling Loss 

-2. .5 


dB 

Antenna Gam 

13 


dB 

Path Loss 

-168. 5 

-169. 7 

dB 

Polarization Loss 

-3 

-4 

dB 

Satellite Receiver 

Antenna Gain 

8 


dB 

Coupling Losses 

-1. 8 


dB 

Link Loss 

’154. 8 

-157. 0 

dB 

Antenna Temperature 

800 


°K 

Receiver Noise Temperature 

400 

1200 

°K 

Effective Receive r Temperature 1200 

2000 

° K 

Receiver Noise Power Density 

-197. 8 

-195. 5 

dBw /Hz 

Receiver Noise Power Density /Link 

-43 

-38. 5 

dBw /Hz 

Loss 

5 X I0" 5 

1.4 X 10" 4 

watts /Hz 


Figure 4-7. Noise Power Density /Link Loss (N s /Lg S )~ GCC to Satellite 

(Frequency = 149. 2 MHz) 


4.3. 2.3 Link from GCC to ATS-3 Satellite and -Return 

The effective noise power ratio per ranging tone is* 


[N" 

— A 

N g , N s 1 

[p_ 

- 4 
3 

P sg L gs P g L gs 


Inserting nominal values gives: 



2.9 X 10'^ per Hz, or -45.4 dB/Hz. 


Worst-case values give: 


— = 2, 2 X 10“ 4 per Hz, or —36. 5 dB/Hz. 

Ms 
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4.3. 2.4 Link from GCC to ATS-1 Satellite and Return 


This link will have slightly less margin than the link described in 
Paragraph 4 . 3 . 2. 3 . 

4.3.3 Channel Allocations 

The factors to be considered in determining the channel 
allocations for this experiment are 

(1) An interfering signal has been noted m the ATS-3 
satellite output at approximately channel 28 frequency 

(2) The GCC self-check tones are transmitted in channel 21 

(3) The ATS-3 pilot tone will be transmitted in channel 12 

(4) It is preferable not to allocate the terminal transponders 
any existing OPLE PEP channels. 

Hence, suitable transmit channels for the terminal transponders 
are OPLE channels 16 and 34. 

4.3.4 Radio Frequency Interference 

A factor which must be considered in the VHF ranging 
experiment is that of radio frequency interference (RFI). Presumably, 
the system will not be purposely jammed, however, the VHF band 
employed by the experiment will be shared with other users. Because 
of the very narrow post detector bandwidth at the Ground Control Center, 
it is unlikely that RFI will appreciably affect the phase measuring cir- 
cuitry. On the other hand, a strong interference signal at or near the 
pilot tone or A/R tone frequency could result in loss of phase-lock, 
either at the aircraft terminal or the Ground Control Center. However, 
such an interference signal can be detected at the GCC by means of the 
satellite monitoring circuitry. 

Because of the heavy use of the VHF band of interest by 
commercial air traffic, it is highly probable that interference problems 
will be encountered during the performance of the VHF Satellite Navi- 
gation Experiment. In this event, two alternatives exist. First, the 
GCC and terminal transponders could be modified to permit realloca- 
tion of channels and thus place the interference signal in an unused 
channel. Second, the experiment could be rescheduled for periods of 
minimum interference. 

4.4 Ionospheric Compensation 

The VHF range estimates used in this experiment contain 
significant bias errors due to the ionosphere. In the following sections. 
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the magnitude of this bias error and a method for introducing a bias 
correction term are discussed. 

4.4.1 The Need for Ionospheric Compensation 

The estimate of the range between the aircraft and a satellite 
contains a bias error due to the ionosphere and the troposphere. The 
magnitude of these bias terms is shown in Figures 4-8 (Reference 2) 
and 4-9. The tropospheric bias is independent of frequency but is a 
function of the temperature of the troposphere, the total pressure, and 
the partial pressure of the water vapor. Even when these quantities 
vary widely the tropospheric bias does not depart significantly from the 
curve in Figure 4-8. Hence, if the elevation angle of the aircraft to the 
satellite is known, or even roughly known, a correction factor can he 
obtained from Figure 4-8 to remove the tropospheric bias. This correc- 
tion is of a small magnitude (35 meters and 3 m3 53 4s for 3® and 37° 
elevation angles, respectively) and can be made with a high degree of 
accuracy. 

On the other hand, the ionospheric bias is shown in Figure 4-9 
to be of significant magnitude and therefore deserves closer attention. 
For example, daytime values of integrated electron density typically 
vary between 2 X 10^7 and 6 or 8 X 10^ electrons /rr?; hence, the iono- 
spheric bias is as large as 1.5 Km to 5 Km for 3° elevation angles. 

When the elevation angle increases to 37°, the bias is still considerable 
(0.6 Km to 2,5 Km). During periods of increased sunspot activity, the 
integrated electron density sometimes exceeds 1 X 10*® electrons /ir?. 

In this worst-case condition, the bias becomes as large as 8 Km. 
Needless to say, some method is needed to remove this bias from the 
range estimate. Various possibilities of doing this are discussed in the 
next paragraph. 

4.4.2 Possible Methods of Ionospheric Compensation 

The methods for measuring or estimating the ionospheric bias 
may be classified as direct and indirect measurements. The former 
measures the ionospheric bias by taking the difference between biased 
and unbiased range estimates while the latter measure s the integrated' 
electron density to determine the ionospheric bias. 

4,4.2, 1 Direct Ionospheric Bias Measurements 

This measurement is made by taking the difference between 
biased and unbiased range estimates of the distance between the Ground 
Control Center and the ATS-3 satellite. The difference between the two 
measurements is attributed to ionospheric bias. Of course, careful con- 
sideration must be given to equipment calibration and accuracy to insure 
that other effects do not add to the ionospheric bias. 
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TROPOSPHERIC BIAS (KM) 



Figure 4-8. Tropospheric Bias Versus Elevation AngLe 
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MAGNITUDE OF IONOSPHERIC BIAS ON ONE-WAY PATH (KM) 



Figure 4-9. Magnitude of Ionospheric Bias Term as a Function 
of Integrated Electron Density 
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The first of two direct bias measurements uses range 
estimates made at YHF and SHF. Because the ionospheric bias decreases 
as the square of the frequency, the bias remaining in a range estimate 
made at 4, 5 GHz wouLd be 1 /900^ as large as the bias in the VHF mea- 
surement at 149.22 MHz. This reduces the bias to a few meters which 
can be neglected when compared to bias of several kilometers in the 
VHF range estimate. The advantage of this method is that it gives an 
accurate result, the accuracy being related to the SHF frequency. It 
suffers the disadvantage of requiring a measurement by other equipment 
(l. e. , SHF equipment) which is not an integral part of the Ground 
Control Center, 

The second method for direct bias measurements uses VHF 
range estimates of the ground station/satellite distance made during 
the night and at the moment of interest. Between midnight and two a. m. 
local time, the ionospheric bias reaches its minimum because of the 
nighttime reduction m ionization of the ionosphere. Hence, a VHF range 
estimate made during these hours will have a minimum bias compared 
to those made of other hours of the day. Typical night values of integrated 
electron density are 3 to 6 X 10^6 electrons/m^ which correspond to 
biases of 100 to 350 meters (for elevation angLes from 3° to 37-). 
Assuming that the true range to the satellite does not vary over several 
nights (a good assumption), the minimum of several nighttime range 
estimates could be used as the reference. This reference is then sub- 
tracted from the satellite range estimate made at the moment of interest. 
The difference, although less accurate than the first method, is an indi- 
cation of the ionospheric bias along the ground station/satellite path. 
However, the easy availability of the reference measurement at night 
makes this approach very attractive. Therefore, because of its rea- 
sonable accuracy and easy implementation, this second method of mea- 
suring the ionospheric bias is the approach to be taken m this 
experiment. 

4.4. 2.2 Indirect Ionospheric Bias Measurements (Faraday Rotation) 

It is well known that the angle through which a polarized wave 
is rotated in traversing the satellite /ground station path is directly pro- 
portional to the integrated electron density of the ionosphere. Hence, by 
measuring the angle of polarization of the received wave at the ground 
station, the ionospheric bias can be determined. However, this method 
has the disadvantage of requiring the use of supplementary antennas and 
receivers to make the Faraday rotation measurement. The need for 
additional equipment makes this approach unattractive from the view- 
point of the present experiment. 
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4,4.3 Correlation Between Measured Ionospheric Bias and ’ 

Predicted Ionospheric Bias \ 

The ionospheric bias is measured for the ground station/ 
satellite link m all the cases discussed above. However, the range 
estimate of the satellite /aircraft distance requires that the iono- 
spheric bias be known for that link also (see Appendix B) . Therefore, 
the important question is how does the ionospheric bias at one location 
relate to that at another location. There is extensive research being 
conducted on this problem by Dr. Aldo da Rosa of Stanford University. 
The desired result of his work is a quantitative relationship between 
the integrated electron density at two different locations. The important 
parameters in this relation are the time of day, latitude, temperature 
profile of the ionosphere, and transition height between oxygen and 
hydrogen molecules m the atmosphere. Dr. da Rosa has made daily 
measurements of the integrated electron density at points separated by 
a few hundred kilometers and at points separated by several thousand 
kilometers. Insufficient data is available to really test his model 
thoroughly. 

In the meantime, some preliminary observations can be made 
regarding the variation of the integrated electron density as a function 
of location. The principle variation is a latitude effect. This latitude 
effect is closely associated with the variation m the F^ layer electron 
concentration as a function of latitude. J. O. Thomas (Reference 3) has 
studied this variation with latitude and has plotted the variation as a 
function of magnetic dip angle for several months through the fall and 
winter of 1957 (see Figure 4-10). Figure 4-11 shows how the median 
of this latitude effect varies over a day in December 1957 Because 
these curves are plotted as a function of magnetic dip angle, one must 
know the magnetic dip angles for the locations of interest. 

Figures 4-12 (Reference 4) and 4-13 (Reference 5) show the 
diurnal variations m integrated electron density at Hawaii /Stanford 
and Stanford/Ely, Nevada. The Hawaii /Stanford data was taken from 
Syncom III in 1964 while the Stanford/Ely data is more recent and was 
obtained from the ATS-1 satellite. The magnetic dip angles for Stanford 
and Ely are 63° N and 64° N, respectively. From Figure 4-10, the 
latitude effect between these two points is seen to be negligible as borne 
out by the close correspondence of the curves in Figure 4-13. Hence, 
for locations which are within a few hundred kilometers of one another, 
the integrated electron densities are essentially identical. Dr. da Rosa 
has much more Stanford /Ely data to substantiate this result. 

However, Stanford and Hawaii have magnetic dip angles of 
63° and 40°, respectively. This results in a significant difference in 
the levels of integrated electron densities over the -measured paths 
because of the latitude difference between the two locations. 
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Figure 4-10. The Variation with Magnetic Dip of the Monthly Median 
Local Noon Values of N m Fg for Different Months 



Figure 4-11. The Variation with Magnetic Dip of the Monthly Median 
Values of N m Fg at Different Hours of the Day in December 1957 
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Figure 4-12. Comparison of Columnar Electron Content Versus Time 
Curves Obtained from Syncom III at Hawaii and Stanford (The Stanford' 
Curve is Shifted to the Right to Correspond to Local Time at Hawaii) 
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COLUMNAR ELECTRON CONTENT < t 0 1 7 ELECTRONS/M 2 ) 



Figure 4-13. Comparison of Columnar Electron Content 
Versus Time Curves Obtained from ATS-1 at Ely and 
Stanford. (The Ely curve is shifted approximately 
one-half hour to the right to correspond to local time at Stanford. ) 
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The 1957 data in Figure 4-10 indicate that the Hawaii levels should be 

1.5 8 and 1.3 8 times the Stanford data in September and October, 
respectively. These figures should be closer to 2 and 2. 5 to more nearly 
relate the data in Figure 4-12. At the moment, the latitude effect curves 
for the appropriate days in September and October 1964 are not avail- 
able. However, Dr. da Rosa expects to have prediction formulas avail- 
able soon that will relate even more closely the Stanford/Hawaii 
observations (Reference 6). Therefore, for this experiment the iono- 
spheric bias as measured at the Ground Control Center wlII be used 
together with appropriate elevation angle and latitude effects to correct 
the aircraft/satellite range estimates. 

4. 5 Velocity Considerations 

The major effect of aircraft velocity is to shift the carrier 
tone frequency received by the GCC. Because this carrier tone is 
tracked by phase-lock loops in the GCC, the doppler shift frequency 
must be limited to the error frequency tracking range of the phase -lock 
loops. This has been demonstrated to be greater than ±300 Hz. There- 
fore, the maximum allowable aircraft velocity in this experiment is 
1170 knots, a value which produces a frequency shift of 3 00 Hz on the 
transmitted 149.2-MHz carrier. This value is set assuming that all 
other carrier frequency errors are compensated for in the GCC through 
calibration settings of the two VHF receiver synthesizers, one for each 
satellite. The allowable doppler shift frequency range could be expanded 
by redesign of the phase-lock loops m the GCC. 
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5.0 


POSITION LOCATION 


There exist several coordinate systems, each with their own 
set of equations describing the satellites and aircraft geometry, which 
can be used in discussing the position location of the aircraft terminal. 
However, there is one set of coordinates with its associated equations 
that is most convenient in the sense that the equations can be solved 
in a straightforward manner. This is the set of orthogonal earth 
centered coordinates together with radius vectors defined in terms of 
direction cosines. It is the purpose of this section to discuss the static 
geometry of this position location experiment. 

The geometrical relation between the satellites and the air- 
craft terminal is shown in Figure 5-1. The Ground Control Center is not 
shown because it does not enter directly into the position location 
geometry. However, it does enter into the earlier stage of measuring 
the ranges Rj and Rg . See Figure B-l and Appendix B for the way m 
which Rj and Rg are estimated from the sidetone phase measurements. 

In this experiment, the satellite locations are assumed to be 
known. A precise statement of the problem is 



GIVEN 

i-H 

O 

Oh 

= RQj^bj = vector from center of the Earth to 
the ATS-1 satellite 



R 03 

= R()3 bg = vector from center of the Earth to 
the ATS-3 satellite 

where 


T 




nT 

= (To j ^ » ^12’ R 13) * s a- vector whose elements are 


the direction cosines for the ATS-1 satellite, 
and bg is the similar vector for the ATS-3 
satellite . 

MEASURED R^ = range between the ATS-1 satellite and the 

aircraft terminal 

Rg = range between the ATS-3 satellite and the 
aircraft terminal 

R^ = range between Earth's center and the air- 
craft terminal. 

ESTIMATE: = R^ = the position vector of the aircraft 

te rminal 
T 

b = b^^j "b a g ) = the direction cosine vector 

for the aircraft terminal. 

i 
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Figure 5-1. Geometrical Relation Between the Earth, Satellites, 

and Aircraft Terminal 


The unknown is b a , a three -component vector. To determine these 
components, three independent equations are required. The law of 
cosines yields two equations which are 


E 1 - R 01 + R h _ 2R 01 R h£a' £l I 5 - 1 ) 

r 3 2 = r 03 2 + R h + R h " 2R 03 R hb^‘b 3 - (5-2) 


The third equation is the auxiliary equation for the direction cosines: 

2 2 2 
b al + b a2 + b a3 = 1 * 


Equations (5-1) and (5-2) yield: 


R 01 2 + E h 2 - R 1 2 


b a b l - b ll b al + b 12 b a2 + b 13 b a3 = 2R 01 R h - Dj (5-4) 

. T V 2 + R h - R 3 2 

— a —3 = b 31 b al + b 32 b a2 + b 33 b a3 = 2R n ,R, “ D 3‘ (5 ~ 5) 

uo h 
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The three independent equations to be solved for the direction cosines 
of the aircraft terminal are (5-3), (5-4) and (5-5). However, the non- 
linear equation in (5-3) complicates the solution somewhat. 

The present satellite positions are approximately 


ATS-1 155° W longitude, 0° latitude 

ATS -3 75° W longitude; 0° latitude. 

The satellites are in an orbital plane that is inclined 0.3° to 
the equatorial plane. The uncertainty in the positions of the satellites is 
±2 Km with the maximum of the error ellipse perpendicular to the 
orbitaL plane. As pictured m Figure 5-1, longitude is measured from 
the X 2 axis toward the Xj axis. The direction cosines are measured 
from each of the three orthogonal axes to the vector in question. For 
example, bgq is the cosine of the angle between the X^ axis and the posi- 
tion vector Fj 1 . Using this notation, vectors bj and b^ are: 



0.423 


cos 65' 


cos 1 5° 1=1-0. 905 


cos 90°— 

cos 1 5 0 
cos 75° 
cos 90 °-</> 



(5-6) 


(5-7) 


and (j ) 3 are small angles which account for the fact that the orbital 
plane of the satellites is tilted a maximum of (5. 24±0. 0475)(10 -8 ) radians 
with respect to the Earth’s equatorial plane. The error term represents 
the uncertainty m satellite position. 

Because the angles 0 ^ and are small, they can be neglected 
on a first-cut at the solution. The steps involved in determining the 
position of the aircraft terminal from Equations (5-3) to (5-5) are 


STEP 1 : Assume b^ = ^33 = 0 

Equations (5-4) and (5-5) become: 


b H b al + b l2 b a2 = B 1 (R 1» R h» R 0 1 ) ( 5 ‘ 8 ) 

b 31 b al + b 32 b a2 = D 3 ( R 3 ’ R h> R 03 ) 
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from which 


STEP 2- 


STEP 3: 


b 32 D l “ b 12 D 3 

b al = Z 

b U D 3- b 31 D l 
b a2 = A 

where 


A - h n b 32 - b 31 b 12 . 
Also, from Equation (5-3): 


I 2 2 

b a3 = ± v 1 ~ b al “ b a2 • 


(5-10) 

(5-11) 


(5-12) 


The Transformations to Latitude and Longitude 


Latitude = 90° - cos“^b a3 

Longitude = tan“l 
Correction for Satellite Position 

Solving Equation (5-4) for b a ^ when hj 3 ^ 0 yields: 




(5-13) 

(5-14) 


(5-15) 


where b a j and b a3 were previously defined in Equations (5-10) 
and (5-12). Likewise, solving Equation (5-5) for b a2 when b 3 3 
^ 0 yields: 




~ b 31 b al 
b 3 2 


b 23 b a3 
b 3 2 



(5-16) 
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The effect of considering the small angles (j) ^ and is to 
change the direction cosines. The third direction cosine is 
likewise affected: 



(5-17) 


The new latitude and longitude which result from the correc- 
tion are obtained by using the new direction cosines in 
Equations (5-13) and (5-14). 


STEP 4 : Correction for Range Errors 

Range errors ARj and ARj result in changes in and 


Ri A R i R 3 AR 3 

AD, = -= , and AD- =— — . 

1 R 01 R h 3 R 03 Rr 

These errors result in changes in the direction cosines 
as given by: 

b 3 2 ADj - b i2 AD 3 
al ~ al A 


(5-18) 


(5-19) 



= b 


a2 


+ • 


’ll 


ADg — bg J AD 


(5-20) 


This carries over to a change in b a g and then into the final 
latitude and longitude of the aircraft terminal. These errors 
and their effects are examined m greater detail in Section 6.0. 
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ERROR ANALYSIS 


PRECEDING PAGfc Blank 


not filmed. 


6 . 0 


Two types of errors are of primary importance m an error 
analysis of this experiment. Namely, they are the measurement errors 
and the final position location errors. The measurement errors refer 
to errors or uncertainties in the measured quantities which are used in 
the position location formulas of Section 5.0. These include the two 
airc raft/ sate Hite ranges, Rj and Rg, the aircraft /Earth center range, 
R^, and the satellite position. Sections 6.1, 6.2, and 6.3 discuss the 
sources of error and the magnitude of the error variance for each of 
the above cases. 

The position location accuracy covered in Section 6.4 relates 
the errors m Rj, R 3 , R^, and satellite position to errors in the posi- 
tion location. The position location error is the important final con- 
sideration, but it is so closely related to the range errors and satellite 
position uncertainty that an error analysis is not complete without con- 
sidering both aspects. 


6 . 1 YHF Range Measurement Errors 

In determining the aircraft's position, one must know the 
aircraft to ATS-1 and aircraft to ATS-3 range in addition to the air- 
craft altitude. These ranges are referred to as and R 3 , respectively 
(see Appendix B, Figure B-l). Because these ranges cannot be mea- 
sured directly, they are estimated by appropriately combining several 
measurable ranges. The relationships are: 


R 1 R g3alg ~ ly/ 2 ( R g3a3g + R glg* 
R 3 = 1 / 2 < R g3a3g ~ R g3g)- 


( 6 . 1 - 1 ) 
( 6 . 1 - 2 ) 


In terms of the sidetone phase delay measurements, the range estimates 
are: 


and 


R 1 = 


Af 21 


09 1 — 1/2001 

c 1 g3 a 1 ~ L 1 

'c^ s txop 


g3alg ' ^ J 'g3a3g 
b 


- 1/2 021 


g*g 


- csc^As frnnC - P 2 J «c K c 


(6.1-3) 


R, = 

3 2 Af 2i 


® 21 g3a3g - ®21 


g3g 


cs Ci // B As tropB -p 2 Q: B K B - 


(6.1-4) 



The mean error in R} and R 3 is non-zero, resulting from 
non-zero means in the tropospheric and ionospheric bias correction 
terms. The mean square error in Rj and R 3 is the sum of the mean 
square errors due to calibration, phase noise, ionospheric bias, and 
multipath under the assumption that each of these errors is independent. 
Our main interest is in the error variance, which is the difference 
between the mean square error and the square of the mean error. 
Because of the independence of the various errors, the error variances 
from each source add to form the total error variance for Rj or R 3 . 

In Sections 6 . 1. 1 through 6 . 1.4, each of the error sources are 
discussed. There it will be seen that the ionospheric error variance 
dominates the other error variances. That is: 

var (ARi) = vax (calibration) + var (phase noise) + 
var (ionosphere) + var (multipath) 

= var (ionosphere) 
and var (AR 3 ) var (ionosphere). 

Because the error variance due to the ionosphere is the same for either 
or R 3 , var (ARi) » var (AR 3 ). This result will simplify later portions 
of the error analysis. 

6.1.1 Calibration Error 

To compute a distance value in sidetone ranging, the portion of 
sidetone phase shift resulting from the delay of propagation through the 
transmission media must be separated from the phase shifts resulting 
from other sources. In this experiment, the effect of other sources of 
phase shift is removed by an off-line calibration process. Residual 
phase shift values not removed by calibration are defined as the system 
calibration error. If daily calibration runs are made, the expected one- 
sigma value of the system calibration error in the VHF range measure- 
ment is 0.3 Km. Major sources of phase shifts which are removed by 
calibration are the tone filters and phase detectors of the GCC, the 
filters in the satellite transponders, and the filters in the aircraft 
terminal. 

The calibration procedure for this experiment is to connect 
the GCC directly to the ATS-3 VHF transponder simulator and record 
the resulting zero-range phase measurements. These phase measure- 
ments are the calibration values of phase shift which are subtracted 
from succeeding phase measurements in the off-line processing. 
Residual phase shift sources not removed by this calibration are the 
linearity of the phase detector and analog-to-digital converter m the 
GCC, quantizing error in the digital processing, and the sidetone phase 
shift difference between the ATS-3 satellite transponder and the satellite 
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transponder simulator. A similar calibration procedure with the ATS-1 
VHF transponder simulator is required if the phase shifts of the side- 
tones through the ATS-1 satellite transponder are significantly different 
from the phase shifts through the ATS-3 satellite transponder. 

The phase shift through the aircraft terminal is likewise 
calibrated by connecting it to the ATS-3 VHF transponder simulator. 

With this connection, the zero -range values in the link from the GCC 
to satellite to aircraft terminal and return are measured and used for 
calibration. 

6.1.2 Phase Noise 

Phase noise is that random noise component which is included 
with the side tone signals at the YHF receiver output of the GCC. The 
resulting signal-to-noise ratio of the received sidetones is related to 
an RMS range error as shown in Figure 4-2. At the worst-case value 
of signal-to-noise ratio predicted for this experiment (3 0 dB within the 
1-Hz processor bandwidth), the RMS value of YHF range measurement 
error due to phase noise is 0.4 Km, 

Another consideration of phase noise is its effect on resolving 
the VHF range ambiguity. The probability of incorrect ambiguity resolu- 
tion as a function of signal-to-noise ratio is shown in Figure 4-3. At the 
predicted worst-case signal-to-noise ratio of 30 dB, the probability 
of incorrect ambiguity resolution is less than 10-7, therefore becoming 
a negligible consideration. 

6.1.3 Ionospheric Error Variance 

As mentioned in Section 4.4, a correction term will he 
introduced into the range estimate to account for the ionospheric bias. 
However, there will always be some difference, (3AK, between the 
actual and estimated values of ionospheric bias (Note; Pg = (3g b/2 
(3 = Pg b/2a c ). This difference, or error, contributes to the total range 

error. In fact, its variance adds to the error variances from the other 
sources because of the independence of the errors. 

The magnitude of var (ionosphere) depends very intimately 
upon the manner in which the ionospheric bias is measured at the 
Ground Control Center and extrapolated to other locations. The tech- 
nique to be used m this experiment is discussed in Section 4.4.2 and 
4.4,3, With the present understanding of the ionosphere, a reasonable 
estimate for the error variance in Rj and R-^ is 1 Km^. Of course, this 
error varies with distance between the GCC and the ionospheric links. 
Until more is known about the ionosphere, this error variance will be 
the dominant error. It is hoped that more can he learned m order to 
reduce this error during the experiment. 
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Figure 6-1. Bottom-to-Top Rejection Ratio of DMC33-2 Antenna 

(Horizontal Polarization) 


6.1.4 Multipath Effect 

The multipath effect is a potential source of range error 
at VHF because it is not practical to suppress the multipath signal 
significantly in the aircraft terminal antenna. This is especially true 
at low elevation angles. Clearly, a parameter directly related to the 
magnitude of the effect is the relative attenuation (a) of the multipath 
signal with respect to the direct signal. The magnitude of a is deter- 
mined by the reflection coefficient of the multipath signal at the ground, 
and by the aircraft antenna pattern bottom-to-top rejection ratio 
because the satellite signal arrives above the horizon and the multipath 
signal arrives from below. A plot of this rejection ratio for the DMC33-2 
Antenna (proposed for this experiment) when actually mounted on an air- 
craft is shown in Figures 6-1 and 6-2 for horizontal and vertical polari- 
zation, respectively. 

At low angles of elevation, where the antenna rejection is low, 
the contribution to the attenuation by the multipath reflector is almost 
negligible over the sea, however, it may be significant for vertical 
polarization over land. Hence, a reasonable estimate of a for both 
polarizations over land is given in Figure 6-1. 
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Figure 6-2. Bottom-to-Top Rejection Ratio of DMC33-2 Antenna 

Vertical Polarization) 


The multipath effect, as it affects ranging, need only be con- 
sidered for range measurements between the satellite and the aircraft. 
The ground station antenna will be sufficiently low to the ground such 
that the differential range between the GCC and the satellite is incon- 
sequential. However, the multipath ray may still cause cancellation at 
the GCC antenna, especially at low angles of elevation. 

Consider the case for the satellite transmitting a carrier at 
u Q and one sidetone at (w 0 +Aw^). Then: 

v ss = cosw t + cos (w 0 +AU|)t . 

As received at the aircraft, the composite direct and multipath signals 
will be: 


v 


ra 


cos(u Q t + 0 od ) + «cos (w Q t + 0 orn ) + cos 


(w 0 + Au^t 
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where a is the relative attenuation of the multipath signal, the subscript 
d refers to the direct path, the subscript m refers to the multipath, and 
the phase angles are related to the ranges and frequencies by: 


R 

e od “ “ u o 

0 ld = c < u o + Au l) 

R 4 AR 

^om - c u o 

R 4 AR , , „ , 

9 lm = c . (U ° +A “1 ) 


Because each tone is finally limited and the phase difference measured, 
the carrier and sidetone should each be expressed as the resultant of 
the direct and multipath signals. Hence 


v 


ra 


A cos (w Q t + 0 od 4 A0 O ) 4 B cos 


(u 0 4 AwjJt 4 0^ d 4 A0j . 


The phase shifts shown by A0 are the errors introduced in each case by 
the multipath signal and are given by: 


a s in(0 orn — d od ) qsin A ru Q 

I4acos(0 om -0 od ) ^acosAruiQ 


asin(0 Lm - 0 Id ) asinATtw^Awj ) 

tan A0^ l4acos{0^ m — 0 ld ) Ibcecos A t(w o 4A(0]^) 

where At is the multipath additional time delay and is AR/c, where AR 
is the multipath range difference. 

In retransmission by the aircraft terminal to the satellite, the 
same multipath model can Be assumed to exist because the time delay is 
very small. Hence, the attenuation constant (a) will be unaltered. Thus, 
the phase difference between the sidetone and carrier on the two-way path 
will be just twice the one-way value, and' 


SlO = 2014+ 2A01 - 2« od - 2A9 0 


2RAW2 

c 


4 2 tan*" 1 


asinA^u^Au 

14Q!cosAt{io o 4Au 1 ) 


— 2 tan -1- 


ffSinATU 0 

l4acosATW o 


The difference in the two arctan terms is the error due to the multipath 
effect. It will be noted that even for ct as large as 1, each arctan term 
will only have a maximum value of 90°, and they are nearly equal. 
Hence, the difference will he small. Another significant fact is that 

is integrated at the GCC over a period of 1 second, and because At 
varies in that time, the error will be further attenuated. 
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To estimate the possible error, therefore, a distribution was 
assumed for At; it was assumed to be Gaussian distributed 'about a 
mean of 1 0 microseconds. The variance in At was taken as 2 X 10“® 
seconds and 2 X 10"^ seconds, corresponding to a differential path Length 
change of 6 meters and 60 meters, respectively. The mean and variance 
of the error term were calculated for f Q = 150 MHz and Afj = 1 kHz, In 
each case, the mean was very nearly zero and the variance was a func- 
tion of a. The results are plotted in Figure 6-3. It is seen that at a = 0.5 
the variance of the error is only 1.5° for a mean delay of 1 0 micro- 
seconds, and this corresponds to a range error of approximately 1 Km 
without smoothing. 

The variance assumed for the differential path length change 
(and hence At) was estimated by observing that the multipath length will 
change in the period of the measurement of 9 -^q (1 second) due to the 
motion of the reflector point across an uneven Earth profile. This is 
especially so at low angles of elevation where a is high. A mean differ- 
ential path length change of 6 meters was considered representative of 
smooth Earth conditions. Because of the relatively rapid fluctuations in 
the reflector profile at aircraft velocities, the small phase error will 
fluctuate rapidly and hence will be further attenuated by the integration 
used m the estimation of Siq* Hence, it is concluded that for elevation 
angles to the satellite m excess of 5°, the multipath error will be 
negligible. At lower angles, the error may become significant, the final 
effect being determined by the spectrum of the error term as well as the 
factors considered here . 


6.2 Altitude Measurement Error - 

The measurement of R^, the distance from the Earth's center 
to the aircraft terminal, is the sum of the Earth's radius and the air- 
craft altitude. Each of these terms include errors as discussed in the 
Appendixes. The variance m the aircraft altitude measurement is 
approximately 160 m2, and the variance in the estimate of the Earth's 
radius is about 100 m.2. The total variance in is thus 260 which 

is a negligible error in comparison with the errors in R^ and Rg, thus 
it will be neglected in the final position location error analysis. 

6 . 3 Satellite Position Error 

Although the positions of the satellites are assumed to be known 
in this position location experiment, there is some uncertainty in the 
precise coordinates of the satellites. The orbital plane of the ATS-1 and 
ATS-3 satellites are inclined 0.3° with respect to the equatorial plane. 
Furthermore, there is a ±2-Km uncertainty in the precise position loca- 
tion of the satellites. The maximum of this uncertainty is perpendicular 
to the orbital plane. 
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RMS PHASE ERROR DUE TO MULTIPATH (DEGREES) 



Figure 6-3. Multipath RMS Phase Error Without Smoothing 
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6.4 


Position Location Accuracy 

All of the previous discussion in Section 6.0 has been centered 
on the errors in range to the aircraft, the aircraft altitude, or the 
satellite position. However, the final and most important question is 
just how these various errors transform to an error m the position of 
the aircraft. In this section, we consider the effect of the range errors 
and the uncertainty of the satellite position upon the final position 
location accuracy. 


6.4.1 Position Location Errors Due to Range Errors 

In Section b. 0, the equations of the position location system 
were derived. Of great interest at this point is the extent to which errors 
m Ej and R 3 affect the final aircraft position location accuracy. In 
Section 5. 0, the aircraft position was defined by a position vector = 

Rj^ b a - -^h ls distance from the Earth's center to the aircraft and is 

the three -component vector of direction cosines defining the aircraft 
position with respect to the Earth's center. As mentioned in Section 6.2, 
the error in is negligible and hence is not considered here. However, 
the errors in Rj and R 3 do deserve further consideration. 


From Section 5.0, the changes m the first two direction cosines 
due to range errors are (from Equations 5. 19 and 5.20): 


where 


db 


b 3 2 0D 1 “ b 12 9 - d 3 


al 


Ob 


b l 1 “ b 3 1 1 


3D 


a2 A 

-Hi -R3 

1 = 9R 1 and dD* = — — 9R 3 . 

"OPh R 03 R h 


(6.4. 1-1) 

(6.4. 1-2) 
(6.4. 1-3) 


It also follows from the algebra of small numbers that the change in 
b a 3 due range errors is: 


3b 


a3 


a3 


a 1 


8b al + b a2 8b a 


2 + 2 ( 9b al 2 + 9b a2 2 ) 


.( 6 . 4.1-4) 


Of great interest is the extent to which these changes m direction cosine 
change the latitude and longitude of the aircraft. Again, from the 
algebra of small numbers, the result is: 


9 (Latitude) = 9b a 3 


(6.4. 1-5) 
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3 (Longitude) = 


(6.4 1-6) 


b al 9 b a 2 ~ b a 2 9b al 

b al 2 + b a 2 2+b a l 9b al + b a 2 ab a 2 

It is seen from Equations 6.4. 1-5 and 6.4 1-6 that the change 
m latitude and longitude due to range errors is a function of location. 

To understand the relationship between range errors and latitude- 
longitude changes, Equations 6.4. 1-5 and 6.4. 1-6 should be pLotted for 
the locus of points covered by the two satellites In lieu of this lengthy 
calculation, a value for a typical position has been calculated A 
latitude of 40° N and longitude of 90° W {Central United States) and a 
range error of 1 Km in Rj and R 3 yield the following results- 

9b al = 1 . 705 X 10 " 4 

3b 9 = - 8 . 25 X 10 " 5 

9b a3 = - 1.435 X 10“ 4 

3(Latitude) = — 1.435 X 10“^ radians 

3(Eongitude) = 1.28X10" 4 radians 

S(Latitude) X R = - 0.915 Km 

3{Longitude) X R = 0. 817 Km. 

This result states that a 1-Km error in range transforms to a 0. 915-Km 
change in latitude and a 0. 817-Km change in longitude. 


Position Location Errors Due to Satellite Position Uncertainty 


The extent to which the ±2-Km uncertainty m the satellite 
position affects the aircraft position location can be seen by examining 
the change in the direction cosines of the aircraft position vector. This 
error is most easily introduced into the problem as an uncertainty in 
the direction cosines (b ^ 3 and bg^) for the ATS-1 and ATS-3 satellites. 
Prom Section 5. 0, the change in the direction cosines due to a change 
in the satellite position is given by: 



(6.4. 2-1) 



(6.4. 2-2) 
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The ±2 -Km maximum uncertainty in satellite position implie s that 
8 t>i 3 = 3b 33 = ±4.75 X 10“^ radians. As in Section 6.4.1, it is seen that 
the position location error is a function of position. For the same Central 
United States location (40 °N, 90°W) as was used, m Section 6.4. 1, the 
results are- 


9b 

9b, 


al 


a 2 

9 (Latitude) 
9(Latitude) X R 
9 (Longitude) 
9( Longitude) X R 


= t7. 23 X 10 “ 5 
= =f1. 06 X 10 “ 4 
= ±1. 51 X 10 " 4 
= 1 . 5 1 X 1 0 ~ 4 radians 
= 0. 962 Km 
= 4.24X10 -5 radians 
= 0.271 Km. 


Again, it must be remembered that these are worst-case errors for 
the ±2-Km error. Checking three points across the United States at 
40°N latitude yields the following latitude /longitude worst-case errors 
due to satellite position uncertainty: 


9 (Latitude) X R 
9( Longitude ) X R 


40 °N /East \ 
75°W\Coast/ 

0. 962 Km 

0. 2705 Km 


40°N /CentralV 
90°W'\U.S. / 

0. 718 Km 

0.462 Km 


40°N /West \ 
120°W\Coast/ 

0. 0942 Km 

0. 737 Km 


It will be possible to remove some portion of these worst-case errors 
when more is known about the motion of the satellite perpendicular to 
the orbital planes. 
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7.0 MANAGEMENT 

The Goddard Space Flight Center (GSFC) of the National 
Aeronautics and Space Administration will conduct the scientific investi- 
gation in the VHF Satellite Navigation Experiment and will produce the 
necessary equipment to demonstrate operational feasibility. 

Responsibility for the investigation is assigned to the Data 
Collection Systems Section of the Communications Research Branch of 
the Systems Division, Goddard Space Flight Centex. Mr. A.E. Jones is 
Chief of the Division, Mr. R.H. Pickard is Head of the Branch, and Mr. 
W.I. Gould is Head of the Data Collection Systems Section. 

Contractors will be utilized to develop and fabricate the 
various equipment necessary for the experiment. Administration sup- 
port is being provided through Dr. Michael J. Vaccaro, Assistant 
Director, Office of Administration. 

The Principal Investigator is Mr. Gay E. Hilton, Senior 
Engineer m the Data Collection Systems Section. He is responsible 
for defining the overall goals of the experiment and for assuring that 
the various equipments are technically capable of meeting these goals. 

The co-mve stigator is Mr. Roger HolLenbaugh, Systems 
Engineer in the Data Collection Systems Section. He is responsible for 
the system study contracts and will be the Technical Officer on one or 
more hardware contracts. He will also be responsible for the Data 
Acquisition and Data Reduction efforts. 

Mr. Walter K. Allen of the Data Collection Systems Section is 
the Project Manager. He is responsible for the overall implementation 
of the VHF Satellite Navigation Experiment, including contract manage- 
ment and coordination of the efforts of the various groups participating 
in the experiment. 

The Tracking and Data Systems {T and DS) Directorate is 
responsible for the operation of the Ground Control Center at GSFC. The 
T and DS Directorate will also support the data analysis requirements of 
the experiment. T and DS support for this experiment will be defined 
at a later date. 
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8.0 SCHEDULE 

Figure 8-1 shows the schedule for the YHF Satellite Navigation 
Experiment. Three months' down time (9, 10, 11) is required of the 
OCC to incorporate the necessary modifications and to perform the 
system tests. These three months should correspond to a lull m the 
present OPLE Experiment, starting after August 1968 and ending before 
February 1969. After modification, the OCC will have the capability 
of collecting data for both experiments, however not simultaneously. A 
short down time in the order of a few minutes will be required to switch 
from one experiment to the other. 
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APPENDIX A 


A COMPARISON BETWEEN RANGING TONES 
AND PN (PSEUDO-NOISE) 

CODES FOR COMMERCIAL AIR 
TRAFFIC CONTROL 


A rather complete comparison between ranging tones and 
ranging codes for range measurements has been made by Baghdady 
and Kruse (Reference 1). Their conclusions are. 

‘'First, in situations, m which there is a decisive requirement 
for high resolution (or precision) and minimum acquisi- 
tion time, with a limitation on allowable bandwidth 
occupancy, systems based purely on harmonic signals 
(or sidetones) are potentially superior to systems based 
exclusively on pseudo-random codes. 

“Second, in situations m which there is a decisive requirement 
for resolving ambiguities in excess of a lunar di-stance, 
systems based purely on pseudo- random codes are 
potentially superior to those based on continuous tones. 

“Third, from the viewpoint of all considerations influenced 
decisively by bandwidth occupancy, systems based 
exclusively on continuous tones are potentially preferable 
to those based exclusively on pseudo-random codes. 

“Fourth, m situations in which there is a decisive requirement 
for transmission security and jam resistance, pseudo- 
random codes are clearly more advantageous. 

“Fifth, in all other respects than those just enumerated, the 
two types of systems are of comparable ratings, “ 

Reviewing these points with respect to this YHF Satellite 
Navigation Experiment, we find that the first point favors the harmonic 
signals from the standpoint of high resolution, minimum acquisition 
time, and minimum bandwidth occupancy. The acquisition scheme with 
the tones is just an algorithm which can be performed very quickly by 
a digital computer as opposed to acquisition by the decoding hardware 
for the pseudo- random codes. If many individual aircraft ranges are 
to be determined, as in an operational system, this savings in acquisi- 
tion time provided by the harmonic signals can become quite 
significant, 
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Because the distances involved in aircraft navigation via 
synchronous satellite are much less than lunar distances, the ambiguity 
resolution of harmonic signals is quite adequate. 

As far as bandwidth and resolution is concerned, Baghdady 
points out that the autocorrelation function of the highest frequency 
sinusoid (f^H^) has a width of 1/2 f^ at the base of each lobe. For 
rectangular pulses of width At seconds, the width of the base of the 
autocorrelation function at the zero crossing point is 2 At. Now, con- 
sider rectangular pulses of width At for which the minimum bandwidth 
is 8/At in order that rectanguLar pulses are preserved. Also, consider 
a sinusoid of frequency f^ = 8 /At Hz, the base of its autocorrelation 
function is 1 /2fj^ = At/16. Because the rectangular pulses have an 
autocorrelation function whose width is 2A t, the sinusoid has a factor 
of 2At/(At/16) = 32 improvement in resolution over that of the pseudo- 
random code which fully occupies the 8/At bandwidth. The sinusoid 
occupies a discrete frequency thereby permitting the adjacent frequency 
spectrum to be utilized for other purposes. 

In commercial air traffic control, secure communications 
are unnecessary, hence harmonic signals can be used. Thus, concern- 
ing the five points made by Baghdady with respect to commercial air- 
craft navigation, ranging by harmonic tones is potentially preferable to 
ranging by pseudo- random codes. 
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APPENDIX B 

ATMOSPHERIC BIAS ERRORS IN VHF RANGING 
1.0 INTRODUCTION 

The sidetone ranging technique described herein measures 
the phase deLay between transmitted and received tones (sinusoids). The 
phase delay is directly proportional to the range between the trans- 
mitter and receiver, thus it yields a range estimate. However, the 
troposphere and ionosphere alter the phase delay from its free space 
value thereby introducing bias errors into the range estimate. It is the 
purpose of this appendix to show the relations between the estimated 
range, the phase delay or shift, and the bias errors. Where possible 
the equations are written to describe the sidetone ranging just as it 
will be done by the Ground Control Center in cooperation with the ATS 
satellites . 

A thorough study of electromagnetic wave propagation through 
the atmosphere under all parametric conditions is a very complex task. 
To simplify this task, atmospheric models (References 1,2, and 3) 
which adequately represent average conditions are employed to estimate 
the effect of the Earth's atmosphere upon range measurements. These 
models indicate that the phase delay of a sidetone Is directly related to 
the index of refraction of the propagation medium. 


2.0 PHASE DELAY DUE TO THE ATMOSPHERE (TROPOSPHERE 
AND IONOSPHERE ) 

The phase path length of a VHF radio path between a ground 
station and satellite is given by (Reference 4): 


P(meters) 



(B - 1 ) 


where p is the index of refraction of the medium along the ray path. 
This index differs from unity primarily in the troposphere (0 to 30.5 
Km) and the ionosphere (85 to 1000 Km). Hence, Equation (B-l) can 
be rewritten as; 


(1) dl + 


= / + Ap trop ) ai+f iuai + f (i 

i J iz J u 

/ z riz /*I4 

(1) dl + / A M-trop ~ / A P-ion ^ 
i J \ U3 


(1 - Ap ion ) dl + f (1) dl 


(B-2) 
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where (1, 12) and (13, 14) denote the lower and upper limits of the tropo- 
sphere and ionosphere, respectively. The first integral on the right 
hand side of Equation (B-2) is just the path length, Rg g) while the other 
integrals can be extended to the limits 1 to 2 because Apj- r0 p and A H- lon 
differ from zero only in the troposphere and ionosphere, respectively. 
The subscript gs on Rg S indicates that the range of interest is the 
ground station /satellite range. Equation (B-2) can now be written as: 

A ^TC P d h-f* *^oz dl 2 < B - 3 > 


P = Rgs + CSC l// 


'/ 


^trop 


dh — esc 


Y 


'A p 1 Q n dh ( B - 4 > 


where dl^ = cscf/'dh, and dig = esc aj dh 

and hq = height of ground or aircraft transmitter above sea level 

hg = height of satellite above sea level. 

The angle ^ is the ray's elevation angle at ground level whereas aq is 
the ray's elevation angle at the middle of the ionosphere. Because of 
the Earth's curvature, they are not the same, but they can be related by 
(Reference 5); 


a. = cos aq = (1 — 0 . 928 cos^) ^ ^ , 


The phase delay in cycles is given by: 

P f_ 
c 


0 (cycles) = — = — P . 


(B — 5) 


(B- 6 ) 


Substituting (B-4) into (B- 6 ) yields; 

9(cycles) = ^ [R gs + (cs c^) A s trop - (cbcOj) A s lon (B-7) 


where 
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trop 


f h 2 

J ^ ^tro] 


dh 
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The Apple fcon-Hartree equation gives Ap.i on as a power series m w - *. 
Consideration of the first-order term is usually adequate, and hence 
As^on becomes: 




(B-9) 


w = radian frequency of signal 

N = electron density of ionosphere (electrons /m^) . 

Ap-^rop 1S independent of frequency but does depend upon the tropo- 
spheric temperature, the total pressure, and the partial pressure of 
the water vapor. 


3.0 RANGE ESTIMATION USING THE PHASE DELAY OF 

THE SIDE TONES 

The ranges of interest in this position location system are 
R]_ and R 3 , the aircraft-to- satellite ranges, as shown in Figure B-l. 

R^ and R 3 cannot be measured directly, but instead are derived from 
a combination of measurements given by: 

R^ = 1/2 (Bg3 a 3g “ Bg3g) (B-10) 

■^1 = -^-gSalg — 1 (Bg 3 a 3 g + Rgi g) (B- 1 1) 

where the subscripts denote the range of interest. Fox example, 

Rg 3 a 3 g is the GCC /ATS-3 /aircraft/ATS-3 /GCC path distance. All of 
these terms refer to a path traversed by the sidetone in making a 
round trip from the GCC to a satellite or aircraft. 

Range Rg 3 g is estimated by measuring the phase delay of a 
signal that is sent to the ATS-3 satellite and returned to the Ground 
Control Center. The round trip phase delay experienced by the refer- 
ence tone (A/R tone) is. 



73 




Figure B-l. Relationship Between Ground Control Center, 

Satellites, and Aircraft 


where 

k’lTJ ~ 149.22 X 2sr X 10^ rad/sec 
W 1D ~ 136.6 X Zjt X 106 rad/sec. 

NOTE; Subscript A refers to path A in Figure B-l, 

The slight frequency offset of the A/R tone from the carrier frequency 
is negligible when calculating or Later, when f 2 jj - f^-g is 

used, an exact difference can be obtained. 

The phase delay of another signal is measured in a similar 
manner. For example, consider the simulated 10.2-kHz Omega signal 
which is placed 426,428057 Hz above the A/R tone in the frequency 
spectrum. Its phase delay is- 
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For ease o£ notation, Let 


Ky 


r h l 


Ndh along path A. 


(B- 14) 


The difference between the phase delays for two distinct frequencies 
is the phase delay for a ranging sidetone at the difference frequency. 
The re suit is 


Af- 


21 


'21 


where 


and 


g3g 


P 2 - 


R g3g + 2csc^ A As tropA -f a A b |3 2 K A 


1 + 1 


L"lU lJ 2U “id“zdJ 


Af 21 - f 2TJ -£ 1U = f 2D - f lD “ 426.428057 Hz. 
The phase delays for the other two ranging sidetone s are 

Af 3i 


whe re 
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31 


'41 


g3g 


g3g 
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•41 


R g 3 g + 2csc ^ A As tropA + a A b ^3 K A 


R g3g + 2csc '/ / A As tropA + °A b h K A 


Af 31 

Af 4i 


P* = 


706. 905276 Hz 
941.428057 Hz 
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1 u 3U u lU 


U 3D W 1D 


(B- 1 5) 


(B-16) 


(B-17) 


(B-18) 


(B- 1 9) 


^4 


1 4 1 


'4U W 1U W 4D W 1D 


w 4U ~ u 3U~ W 1U 


149. 22 X 2it X 10 6 rad/sec 


U 4D ~ U 3D K(J 1D ~ 135. 6 X2ttX 10^ rad /sec . 

The relation between the A/R tone and the sidetones in one of the 
2.5-kHz channels is shown in Figure B-2. 
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Figure B-2. Frequency Spectrum, of One 2. 5 -kHz Channel 


Because the phase measurement of #21 ^ an between zero and 

one cycle, the other two phases, 0-3 1 „ and 9 a .~\ , are combined in an 

g3g g3g 

ambiguity resolution scheme to determine the proper number of whole 
cycles associated with each range measurement. 

The phase delay of the Af 21 ranging sidetone along the path 
B g3a3g 1S: 


g3a3g 


^21 f 

~ ^ R g3a3g + 2csci// A^ s tropA + 2csci^gAs tr0 pg 

+ P 2 b ("a k a + “b k b) 1 • < B - : 


Using (B-10), (B-13), and (B-20) to determine R 3 yields: 


% " 1 / 2 ( R g3a3g “ R g3g) 


A 

R 3 = 


2 % 


( 021 g3a 3g 021 g3g) csc ^B AStro P B 


- U 

“ 2 “bS 


(B-21) 
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Range estimates of R 3 can also be made using the phase shifts 
at sidetones Af 3 j and Af^. The result is the same as that shown in 
Equation (B-21) for sidetone Af 21 except that the subscripts 21 become 
31 or 41 and J32 becomes (33 or ( 34 , respectively. 

Repeating the above procedure, range estimates for Rj can be 
derived based on Equation (B-ll). The range estimate Rj made by Afgi 
sidetone is 

/v Af?i / \ 

Ri = [e ? i , - 1/2 e 71 - 1/2 e ? , 

1 c \ 21 g3alg Zi g3a3g 21 glg/ 

- esc c As trop c - |3 2 |«c K c ‘ - (B-22) 

The range estimates of R^ based upon the other sidetones (Afg ]_ and 
Af 4 i) follow immediately. 

There are methods of combining the three range estimates 
in such a way that the errors in the new estimate is less than the error 
in any individual range estimate. Furthermore, possibilities exist for 
removing the atmospheric bias. These are discussed elsewhere m the 
text. 
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DISTANCE FROM THE AIRCRAFT TO THE 
CENTER OF THE EARTH 


1.0 INTRODUCTION 

In the VHF Satellite Navigation Experiment, ranging measure- 
ments to an aircraft are made from two geostationary satellites. After 
proper range ambiguity resolution, two spheres of position are obtained, 
one centered at each satellite. A third sphere of position can be obtained 
by determining the distance from the aircraft to the center of the Earth. 

The determination of this third sphere of position is dis- 
cussed here. It is assumed that the true altitude above mean sea level, 
h.£, has been measured accurately with a barometric altimeter or by 
other means and that the geodetic latitude of the aircraft is known. 

2.0 DEFINITIONS (Jeffreys, THE EARTH, Cambridge 

University Pres s, London, 1959.) 

Geoid: The locus of points below the Earth's surface equal 
to the measured height. It coincides with the mean sea 
level surface on the oceans. 

Referenced spheroid The ellipsoid of revolution which most 
nearly approximates the mean sea level surface of the 
Earth. It is characterized by the equatorial radius R e 
and the flattening f, where: 

f-i R ° 

' R e - e p 

and Rp is the polar radius. 

Altitude h t : The height above the geoid (or mean sea level 
surface) and is measured perpendicular to the refer- 
ence spheroid. 

Geodetic latitude <£>: The angle between the equatorial plane 

and the normal to the reference spheroid at the station. 
Geodetic latitude is indicated on maps. 

Geocentric latitude 9 (declination) - The angle between the 
equatorial plane and the radius vector to the station. 
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Vehicle position (0, <$>) ; The geocentric and geodetic latitude 
of the vehicle S are the geocentric latitudes of the sub- 
vehicle point, or the point P, defined by erecting a 
perpendicular at the reference spheroid which passes 
through S. 

3.0 DISCUSSION (Refer to Figure C-l) 

The error involved in calculating the airplane's geocentric 
radius (r) from its altitude by use of the equation: 

r = h + R(0) (C-l) 

where R(0) ~ R e (1 — f sin 2 0) (C-2) 

is due to three causes. These are: 

(1) The approximation involved in describing the reference 
spheroid to first order in f; 

(2) The approximation involved in defining the altitude along 
the radius vector rather than the vertical to the refer- 
ence spheroid; 


82 








(3) The difference between the geoid and a reference spheroid 
with a specified flattening. 

The reference spheroid is given by 


R 2 cos 2 0 , R 2 sin 2 0 
— -f 


R, 


R 2 (l-f) 2 


= 1 


from which 


R = R e (1-f) l-(2f-f 2 ) cos20|- 1/2 . 
Equation C-4 can be written for small f: 


R « R c 


1 - f sin 2 0 - 3/8f 2 sin 2 20 + 0{f 3 )] 


The second-order term in Equation C-5 has a, value of: 

Ah R = 3/8£ 2 R e sin 2 20 = 26.9 sm 2 20 [m] . 


(C-3) 


(C-4) 

(C-5) 


(C-6) 


Thus, the altitude error in describing the reference spheroid to first 
order in f has a maximum value of 26.9m at 0 = 45 degrees. 

The geocentric latitude (0) of a point on the spheroid is related 
to its geodetic latitude ((f)) by: 


tan {(f) — Q) 


I dR 
R de 


R 


e 


R 


f sin 2 6 + ^f 2 
4 


sin 49 + 0(f 3 ) 


and to first order m f. 


(C-7) 


tarn// - tan {<f> — 9) = f sm20. 
The angle ip has a maximum value of. 


(C-8) 


tf'max = “ 0 ) ma x ~ °- 192 degrees at 9 = 45 degrees. (C-9) 

One can see from Figure C-l that for 9 - 45 degrees,if(>i/(> rj, and the 
maximum altitude error due to the error in the vertical to the reference 
spheroid is approximately: 


Ahjj = h(l — cos tj) ^h(l — cos ip) 



(C-10) 


•with a maximum value of 5.6 X 10 _ 6 at 9 = 45 degrees. For aircraft 
altitudes, this amounts to fractions of a meter and can safely be 
neglected. 
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The difference in height (AIiq) between the geoid and reference 
spheroid depends on the parameters of the reference spheroid. For 
instance, Kanla in 1961 (Reference 1) quotes differences between— 51m 
and +43m, although more recent geodetic and gravimetric data should be 
utilized. 


For the VHF Satellite Navigation Experiment, it is recom- 
mended that the maximum error in airplane radius calculation be limited 
to ±lOm (33 ft) so that the calculations are significantly more accurate 
than the observed data. The expressions for the airplane radius should 
hence be- 

r = h fc + R e (1-f sin 2 -0 - 3/8 f 2 sin 2 20) 4Ah G (C-ll) 

where Iiq can be tabulated and stored for the Earth as a function of 
latitude and Longitude. It is important that the values of R e , f, and AIiq 
be consistent and that agreement be reached on their values. 

Because in general the geodetic latitude of the aircraft is known 
rather than the declination, an approximate value of 0 can be obtained 
from Equation C-8 by noting that: 

tan ( $ - 0 ) ~ <£ — 0 for small { <j> — 9 ) 

and the error introduced in replacing 6 by <j) on the righthand side of 
Equation C-8 is of order f 2 which yields an error of order f 2 in 
Equation C-ll. Thus, Equation C-8 can he written approximately 

0 = 4> - f sin 2<j). (C-12) 

The error m airplane radius introduced by using <j> instead of the more 
correct Equation C-12 is negligible. 

4.0 SUMMARY 

Given values of hj- and <j>', the airplane distance from the center 
of the Earth can be determined to within ±10m by use of the expressions: 

0 = $ 1 - f sin 2 <j) ] 

r = h t + R e (1-f sin 2 9 — 3/8 f 2 sin 2 20) 4 Ah G . 

It has been assumed that the difference (AIiq) between the geoid and 
reference spheroid is known to about ±5m along the flight path of the 
airplane. A detailed analysis of more -recent data is required to deter- 
mine appropriate values of R e , f, Ah.Q ar id the uncertainty associated 
with these values when the VHF Satellite Navigation Experiment will be 
performed. 
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APPENDIX D 

ALTITUDE MEASUREMENT ERRORS 


The primary instrument for measuring an aircraft's altitude 
is the barometric altimeter. Because, measurement of altitude is part 
of the determination of a third sphere of position in the YHF Satellite 
Navigation System, an evaluation of altitude measurement errors is 
required for the estimation of position location accuracy. The aircraft 
is assumed to he a commercial, high-speed, subsonic turbojet and the 
following material has been obtained from Reference 1. 

(1) The system error in measuring altitude is considered to 

be the statistical sum of the instrument error and 
the static pressure error. 

(2) The instrument error is due to drift, friction, temperature, 

hysteresis, and similar mechanical causes as well 
as the interpretation of the altitude and pressure 
setting scales by the pilot. 

(3) The static pressure error is the difference between 

the free-stream static pressure and the pressure 
indicated by the pitot tube or fuselage vent due to 
compressibility and flow effects. The static 
pressure error can be further sub-divided into a 
fixed error for a given aircraft type (this fixed 
error can be partially compensated for) and a 
random error. 

The precision altimeter instrument error (three-sigma value) 
was estimated in Reference 2 to be 132 ft, as recommended by the 
IATA and ICAO, the specified accuracy of air data compensators is 
usually given as 0. 2 percent of altitude (i.e. , 70 ft at 35, 000 ft). Vari- 
able static pressure errors have been estimated by IATA to be 250 ft „ 
between 30, 000-ft and 50, 000-ft altitudes. These error estimates seem 
to be borne out by flight tests described in Reference 1, and it is assumed 
that they are three -sigma values. 

System errors were obtained by making pressure altitude and 
radio altimeter measurements between two pairs of aircraft on routine 
flights over the North Atlantic at altitudes above 29, 000 ft (Reference 3). 
Data reduction yielded a three-sigma error of 510 ft for a pair of air- 
craft and a calculated error value of 360 ft for a single aircraft. In 
Reference 2, three-sigma system errors of 342 ft were estimated for 
aircraft using uncorrected precision altimeters and 279 ft for aircraft 
using servo-corrected precision altimeters. 
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